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We discuss topological properties of the ground state of spatially homogeneous ensemble of fermions. 
There are several classes of topologically different fermionic vacua; in each case the momentum space 
topology of the vacuum determines the low-energy (infrared) properties of the fermionic energy 
spectrum. Among them there is class of the gapless systems which is characterized by the Fermi- 
hypersurface, which is the topologically stable singularity. This class contains the conventional 
Landau Fermi-liquid and also the non-Landau Luttinger Fermi-liquid. Another important class of 
gapless systems is characterized by the topologically stable point nodes (Fermi points). Superfluid 
3 He-A and electroweak vacuum belong to this universality class. The fermionic quasiparticles (par- 
ticles) in this class are chiral: close to the Fermi points they are left-handed or right-handed massless 
relativistic particles. Since the spectrum becomes relativistic at low energy, the symmetry of the 
system is enhanced in the low-energy edge. The low-energy dynamics acquires local invariance, 
Lorentz invariance and general covariance, which become better and better when the energy de- 
creases. Interaction of the fermions near the Fermi point leads to collective bosonic modes, which 
look like effective gauge and gravitational fields. Since the vacuum of superfluid 3 He-A and elec- 
troweak vacuum are topologically similar, we can use 3 He-A for simulation of many phenomena in 
high energy physics, including axial anomaly. 3 He-A textures induce a nontrivial effective metrics of 
the space, where the free quasiparticles move along geodesies. With 3 He-A one can simulate event 
horizons, Hawking radiation, rotating vacuum, conical space, etc. 
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I. INTRODUCTION 

The physical vacuum is a complicated condensed mat- 
ter ly-Q. At the moment we know only the low-energy 
properties of this substance, i.e. the properties at en- 
ergies much smaller than the Planck energy scale, E <C 
Ep = ^/Tic 5 /G, where G is Newton's constant. We know 
that at the low-energy edge our vacuum has many dif- 
ferent symmetries: U(l) and SU(3) gauge symmetries, 
Lorentz invariance, general coordinate invariance, and 
dicrete CPT symmetry. With increasing energy more 
elements of symmetry are added: SU(2) symmetry of 
weak interactions, probably the GUT symmetry and even 
supersymmetry. When the temperature dicreases such 
symmetries become spontaneously broken. This is the 
traditional point of view, which is supported by the ob- 
servation of the similar symmetry breaking in different 
(many-body) condensed matter systems: superfluids, su- 
perconductors, magnets, liquid and ordinary crystals. 

However there is another and actually opposite point 
of view: all symmetries known in the Universe sponta- 
neously (but without any phase transition) appear at the 
low-energy corner [H . They become more and more pro- 
nounced the lower the energy. The symmetries disappear 
at higher energies when the Planck energy is approached. 
This conjecture is also supported by the condensed mat- 
ter analogy. 

Here we discuss the topological origin of such enhanced 
symmetry in the infrared limit. It is caused by the so- 
called topologically stable Fermi points in the spectrum 
of the fermionic system. We show that if such points 
exist in quantum condensed matter, then at low enough 
temperature this system exhibits fully or partially the 
Lorentz invariance, general covariance, and gauge invari- 
ance. Moreover the collective modes which describe such 
a system at low T all are represented by the chiral rel- 
ativistic fermions, gauge fields, and gravity. All are the 
low-energy phenomena, which are absent at higher ener- 
gies, where the condensed matter has only a very limited 
set of the global symmetries. 

From the second point of view some directions in 
physics look artificial. In particular, since the gravity 
exists only as the infrared phenomenon one should not 
quantize gravity: Only low-energy gravitons can be quan- 
tized [p| . The same concerns the multidimensional string 
theories, which also give rise to gravitation in the infrared 
limit. The Fermi point mechanism does not require a 
high dimensionality for the space-time: The topologically 
stable Fermi point is just a property of the conventional 
3+1 dimensional space-time. 



II. MANIFOLDS OF ZEROES. 

The Fermi point is a particular case of the topologically 
stable manifolds of zeroes in the fermionic spectrum. Let 



us start with the simplest such manifold - the Fermi sur- 
face. 



A. Fermi surface as topological object 

1. Topological stability of Fermi surface 

The Fermi surface appears in the noninteracting Fermi 
gas, where the energy spectrum of fermions is 



B <"» = k 



/'■ 



(1) 



and \i > is the chemical potential. The Fermi surface 
bounds the volume in the momentum space where the en- 
ergy is negative, E(p) < 0, and where the particle states 
are all occupied at T = 0. In this isotropic model the 
Fermi surface is a sphere of radius pf — \/2mJl. 

It is important that the Fermi surface survives even if 
interactions between particles are introduced. Such sta- 
bility of Fermi surface comes from the topological prop- 
erty of the Feynman quantum mechanical propagator - 
the one-particle Green's function 



g = {z-ny 



(2) 



Let us write the propagator for a given momentum p 
and for the imaginary frequency, z = ip . The imagi- 
nary frequency is introduced to avoid the conventional 
singularity of the propagator "on the mass shell" , i.e. at 
z = E(p). For noninteracting particles the propagator 
has the form 



G 



ipo - E(p) 



(3) 



Obviously there is still a singularity: On the 2D hypersur- 
face (po = 0,p — pf) in the 4-dimcnsional space (po,p) 
the propagator is not well defined. This singularity is 
stable, i.e. cannot be eliminated by small perturbations. 
The reason is that the phase $ of the Green's funtion 
G = |G|e J * changes by 2tt around the path embracing 
this 2D hypersurface in the 4D-space (see Fig. pi). The 
phase winding number N = 1 cannot change continu- 
ously, that is why it is robust towards any perturbation. 
Thus the singularity of the Green's funtion on the 2D- 
surface in the momentum space is preserved, even when 
interactions between particles are introduced. 

Exactly the same topological conservation of the wind- 
ing number leads to the stability of the quantized vortex 
in superfluids and superconductors, the only difference 
being that, in the case of vortices, the phase winding oc- 
curs in the real space, instead of the momentum space. 
The complex order parameter \P = l^le 4 * changes by 
2nN around the path embracing vortex line in 3D space 
or vortex sheet in 3+1 space. The connection between 
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Fermi point atp=0. 
The direction of spin makes 
a hedgehog in the momentum space 
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Class (4) 
Topologically unstable 
Fermi line on the p axis 

(c) 



singularity in the Bogoliubov spin 
in Class (4) is removed 
after deformation to Class (2) 

(d) 



FIG. 1. (a) The Fermi surface in the 2+1 momentum space. 
For simplicity the p z coordinate is suppressed so that the 
Fermi surface is the line (po — 0,p — pp) in the 2+1 space 
(po,Px,p y )- This line is a singularity, which is similar to 
a vortex in a real 3D-space: The phase of the propagator 
G = (ipo — {Px + Py — pp)/2m)~ 1 changes by 2-7T around the 
line in the momentum space in the same manner as the phase 
of the order parameter changes by 2-7T around a vortex in 
the real space. In general case the winding number Ni - 
the topological invariant of the Fermi surface - is given by 
Eq.(|4|). (b) Fermi point at p = in the 3D momentum space 
(Px,p y ,Pz)- At this point the particle energy E = cp is zero. 
A right-handed particle is considered with its spin parallel to 
the momentum p, i.e. s(p) = (l/2)p/p. The spin makes a 
hedgehog in the momentum space, which is topologically sta- 
ble. The topological invariant of the Fermi point is N$ — 1, 
where N3 is given by Eq.(|13|). (c) Fermi line - topologically 
unstable manifold of zeroes - is shown in the 3D momentum 
space {Px,p y ,Pz). The (Bogoliubov) spin (arrows) is confined 
into the (p x ,Py) plane and has a singularity on the p z axis, 
(d) This singularity can be removed bya continuous transfor- 
mation. The spin escapes into a third dimension (p z ) and be- 
comes well defined on the p z axis. As a result, the quasiparti- 
cle spectrum becomes fully gapped (the "relativistic" fermion 
acquires the mass). 



Collective modes of Fermi systems 



Oscillations in the shape of Fermi surface 




Oscillations of slopes and shifts of Fermi points 





Mode in which position 
of node oscillates 

is the effective dynamic 
electromagnetic field 



Modes in which 
the slopes are oscillating 



E 2 = cV -» S' k PiPk 



is the effective gravity field 





FIG. 2. Collective modes of fermionic systems with Fermi 
surface and Fermi point. 



The topology in the 4+4 dimensional space describes: 
(i) The momentum space topology in the homogeneous 
system; (ii) The defects of the order parameter in real 
space; and (hi) Topology of the energy spectrum within 
the topological defects || . 

In the more complicated cases, when the Green's func- 
tion is the matrix with spin and band indices, the phase 
of the Green's function becomes meaningless. In this case 
one should use a general analytic expression for the in- 
teger topological invariant which is responsible for the 
stability of the Fermi surface: 



N 1= Tr 



f 2ri,^( Po,P ^ dlS 1 ^°' p ) 



(4) 



Here the integral is taken over an arbitrary contour C in 
the momentum space (p,Po), which encloses the Fermi 
hypersurface; and Tr is the trace over the spin and band 
indices. 



the topology in real space and the topology in momen- 
tum space is, in fact, even deeper (see e.g. Ref. M). If 
the order parameter depends on space-time, the propa- 
gator in semiclassical aproximation depends both on 4- 
momentum and on space-time coordinates G(po, p, £, r). 



2. Landau Fermi liquid 

The topological class of systems with Fermi surface is 
rather broad. In particular it contains the conventional 
Landau Fermi-liquid, in which the propagator preserves 
the pole. Close to the pole the propagator is 



G = 



ipa - v F (p-p F ) 



(5) 



Evidently the residue Z ^ 1 does not change the topolog- 
ical invariant for the propagator, Eq.(|4|), wihch remains 
N\ = 1. This is essential for the Landau theory of an 
interacting Fermi liquid; it confirms the assumption that 
in Fermi liquids the spectrum of quasiparticles at low en- 
ergy is similar to that of particles in a Fermi gas. It is 
also important for the consideration of the bosonic collec- 
tive modes of the Landau Fermi-liquid. The interaction 
between the fermions cannot not change the topology of 
the fermionic spectrum, but it produces the effective field 
acting on a given particle by the other moving particles. 
This effective field cannot destroy the Fermi surface ow- 
ing to its topological stability, but it can locally shift the 
position of the Fermi surface. Therefore a collective mo- 
tion of the particles is seen by an individual quasiparticle 
as a dynamical mode of the Fermi surface. These bosonic 
oscillative modes are known as the different harmonics of 
the zero sound. An example is shown in the upper part 
of Fig§ 

Note that the Fermi hypersurface exists for any spatial 
dimension. In the 2+1 dimension the Fermi hypersurface 
is a line - the vortex line in the 2+1 momentum space as 
is shown in Fig. |lj(a). 



3. Non-Landau Fermi liquids. 

In the 1+1 dimension the Green's function looses its 
pole, but nevertheless the Fermi surface is still there 
|j[10|. In 1+1 dimension the Landau Fermi liquid trans- 
forms to another states within the same topological class. 
Example is provided by the Luttinger liquid. Close to the 
Fermi surface the Green's function for the Luttinger liq- 
uid can be approximated as (see ||ll|fl|l2|) 



(ipo 



G(z,p)~ 

9 — 1 9 9 ~ 1 9 

vip) ~^~ {ipo + vip) 5 (ip - v 2 p) ~^~ {ipo + v 2 p) 5 (6) 



where v% and v 2 correspond to Fermi velocities of spinons 
and holons and p = p — pp . The above equation is not 
exact but reproduces the topology of the Green's func- 
tion in Luttinger Fermi liquid. If g ^ and v\ ^ v 2 , the 
singularity in the (p, z = ipa) momentum space occurs 
on the Fermi surface, i.e. at (pa = 0, p = 0). The topo- 
logical invariant in Eq.(??) remains the same N\ = 1, as 
for the conventional Landau Fermi-liquid. The difference 
from Landau Fermi liquid occurs only at real frequency z: 
The quasiparticle pole is absent and one has the branch 
cut singularities instead of the mass shell, so that the 
quasiparticles are not well defined. The population of 
the particles has no jump on the Fermi surface, but has 
a power-law singularity in the derivative pi . 

Another example of the non-Landau Fermi liquid is 
the Fermi liquid with exponential behavior of the residue 



fll3| . It also has the Fermi surface with the same topo- 
logical invariant, but the singularity at the Fermi surface 
is exponentially weak. 



4- Superconducting transition: from pole to zero 

In this section we do not consider the bifurcations 
which lead to the overall reconstruction of the particle 
spectrum, which occur, say, during the transition to su- 
perconducting state. We assume that the temperature, 
though is small compared to the "Planck" scale Ep = p, 
is still above the superfluid transition: T c < T <C Ep. 
But even in the superconducting region T < T c the wind- 
ing number for the conventional Green's function is pre- 
served: the invariant in Eq.(??) is again Ni — 1. Only 
instead of pole at p = pp in the Green's function of the 
normal state, one has zero at p = pp in the supercon- 
ducting state: 



G=- 2 ZP0 + M(P) A2 ,M( P ) 

pI + m 2 (p) + a 2 Vi ^ 



p 

2m 



p m vp(p-pp) 

(7) 



where A is the gap in the quasiparticle spectrum. This 
demonstrates the stability of the singularity in the 
Green's function, though the Fermi surface as the man- 
ifold of zeroes in the energy spectrum disappears in the 
superconductor. 



B. Fully gapped systems. 

Although the systems we have discussed contain 
fermionic and bosonic quantum fields, this is not the 
relativistic quantum field theory which we need for the 
simulation of quantum vacuum: There is no Lorentz in- 
variance and the oscillations of the Fermi surface do not 
resemble the gauge field even remotely. The situation 
is somewhat better for superfluids and superconductors 
with fully gapped spectra; examples, which provide use- 
ful analogies with Dirac fermions and spontaneously bro- 
ken symmetry in quantum fields, are conventional super- 
conductors Q and superfluid 3 He-B Q. In 3 He-B the 
Hamiltonian of free Bogoliubov quasiparticles is the 4x4 
matrix: 



H = 



(M(p) ca-p 
\ca -p -M(p) 

A 



= T 3 M(p) + cT 1 a-p , (8) 



c= — , U l =E Z = M z (p) + A z , (9) 

PF 

where the Pauli 2x2 matrices a describe the conventional 
spin of fermions and 2x2 matrices r describe the Bogoli- 
ubov isospin in the particle-hole space. The Bogoliubov 
Hamiltonian asymptotically approaches the Dirac one in 



the low p limit. This is however is not the low energy 
limit, which for the typical 3 He-B parameters occurs for 
p close to pf- The low p and low E limits coincide at 
large negative /x, which is not the case in 3 He-B. Nev- 
ertheless the 3 He-B also serves as a model system for 
simulations of many phenomena in particle physics and 
cosmology including experimental verification |lq ] of the 
Kibble mechanism describing formation of cosmic strings 
in the early Universe |17| |. 

Note that even for the fully gapped systems there can 
exist the momentum-space topological invariants, which 
characterize the vacuum states. Typically this occurs 
in 2+1 systems, e.g. in the 2D electron system exhibit- 
ing quantum Hall effe ct llq |; in thin film of 3 He-A (see 
JL9[ and Sec. 9 of Ref. [|20f); and in the 2D superconduc- 
tors with broken time reversal symmetry ||2l| . The quan- 
tum (Lifshitz) transition between the states with differ- 
ent topological invariants occurs through the intermedi- 
ate gapless regime |20(|. 



C. Fermi point 

1. Systems with Fermi points: 3 He- A and Standard Model 

Now we proceed to the topological class which most 
fully exhibits the fundamental properties needed for a 
realization of the relativistic quantum fields, analogous 
to those in particle physics and gravity. 

It is the class of systems, whose representatives are su- 
perfluid 3 He-A and the vacuum of relativistic left-handed 
and right-handed chiral fermions (another example of 
this class in condensed matter has been discussed for gap- 
less semiconductors |E2[). This class is characterized by 
points in the momentum space where the (quasi)particle 
energy is zero. In particle physics the energy spectrum 
E(p) — cp is characteristic of the massless neutrino (or 
any other chiral lepton or quark in the Standard Model) 
with c being the speed of light. The energy of a neu- 
trino is zero at point p = in the 3D momentum space. 
The Hamiltonian for the neutrino - the massless spin-1/2 
particle - is a 2 x 2 matrix 



H = ±ca ■ p 



(10) 



which is expressed in terms of the Pauli spin matrices 
a . The sign + is for a right-handed particle and — for a 
left-handed one: the spin of the particle is oriented along 
or opposite to its momentum, respectively. 

The Bogoliubov matrix for the 3 He-A fermions is 

H = T 3 M(p) + c(a • 4)(ne (1) ■ P - r 2 e^ • p) , (11) 

where e^ and e' 2 ' are real vectors, which in equilibrium 
or ground state are unit orthogonal vectors. The energy 
of the fermions 



El = A/ 2 (p) + c 2 (p x l) 2 , 1 = e< 3 ' = e« x e( 2 ) , (12) 



is zero at 2 points, at p = ±pj?l. 



2. Topological invariant for Fermi point 

Let us show that zeroes in the spectrum of the chi- 
ral fermions and in the spectrum of 3 He-A fermions are 
topologically stable and described by the same topolog- 
ical invariant. Let us again consider the propagator of 
the particle Q — (ipo — 7i)~ l on the imaginary frequency 
axis, z — ipo. One can see that this propagator still has 
a singularity, which is now not on the surface but at the 
points in the 4D momentum space: (po = 0, P = 0) and 
(po = 0, ±pf1) for chiral fermions and 3 He-A fermions 
respectively. An integer topological invariant which sup- 
ports the stability of the point zeroes is expressed in 
terms of the propagator in the following way p0[ : 



N* 



1 



24tt 2 



'/(y A7 



tr / ds~< gd p G- l gd p a- l gd Px g 



(13) 



where a is the 3-dimensional surface embracing the point 
node in the 4-momentum space. 

For the chiral fermions in Eq.([10|) this invariant is 
Ng = ±1, where the sign is determined by the chirality 
of the fermion. The meaning of this topological invariant 
can be easily visualized. Let us consider the behavior of 
the particle spin s(p) as a function of its momentum p in 
the 3D-space p = (p x ,Py-Pz)- For right-handed particles 
s(p) = p/2p, while for left-handed ones s(p) = — p/2p. 
In both cases the spin distribution in the momentum 
space looks like a hedgehog (see Fig. |l|b), whose spines 
are represented by spins: spines point outward for the 
right-handed particle and inward for the left-handed one. 
In the 3D-space the hedgehog is topologically stable. 



3. Spin from isospin. 

For the 3 He-A fermions in Eq.([ll|) A3 = +2, where the 
sign is determined by the position of the node, p = ip^l. 
The topological invariant is twice larger because of the 
double degeneracy of the Fermi point over the conven- 
tional spin of the 3 He atom. For each projection of spin 
one has N3 — +1. Note that the Bogoliubov spin f in 
3 He-A plays the same role as the conventional spin a of 
chiral fermions in Eq. ([IO]) . On the other hand the conven- 
tional spin of the 3 He atom is responsible for the degen- 
eracy, but not for chirality, and thus plays the part of the 
isospin (see also Sec. Ill H). This means that the origin of 
the spin responsible for the chirality of the (quasi)particle 
is fully determined by the matrix structure of the Fermi 



point. In this sense there is no principle difference be- 
tween spin and isospin: changing the matrix structure 
one can convert isospin to spin, while the topological 
charge of the Fermi point remains invariant. 



D. Fermi line 

1. Superconductivity in cuprates. 

The high-temperature superconductors in cuprates 
most probably contain zeroes in their quasiparticle en- 
ergy spectrum. The ARPES experiments p3| show that 
these are four lines in the 3D momentum space where the 
quasiparticle energy is zero or, equivalently, four point 
zeroes in the 2D CuC>2 planes. The high-T superconduc- 
tors thus belong to class of systems with Fermi lines: the 
dimension D of the manifold of zeroes is 1 , which is in- 
termediate between a Fermi surface with D = 2 and a 
Fermi point with D = 0. 

The energy spectrum of quasiparticles near each of the 
4 gap nodes can be written as 



H = TlC X (p x -pW) + T3C V (p y - Pi 



- n (f>h 



(14) 



The "speeds" of light c x and c v are the "fundamental" 
characteristics determined by the microscopic physics of 
the cuprates. This means that the system belongs to 
the same class as 2+1 QFT with massless relativistic 
fcrmions. 



3. Topological instability of Fermi line. 

The lines of zeroes generally have no stability: There 
is no corresponding N2 invariant, which can support the 
topological stability. The singular line in the momentum 
space from which the spines (now the vector f) point 
outward (see Fig. |]c) can be elmininated by the escape 
of the r-vector to a third dimension. This can be ac- 
complished by an operation similar to the folding of an 
umbrella (see Fig. Qd). 

Existence of the nodal lines can be prescribed, how- 
ever, by the symmetry of the ground state. There are 
many nontrivial classes of superconductors, whose sym- 
metry supports the existence of nodal lines in symmetric 
positions in the momentum space Eg ]. The symmetry 
violating perturbations, such as impurities, an external 
magnetic field, etc., destroy the lines of zeroes ||. One 
could expect different types of transformations of these 
lines of zeroes which depend on the perturbation. Im- 
purities, for example, can : (i) produce the gap in the 
fermionic spectrum J27J] (see Fig.[yd), which corresponds 
to appearance of mass for the 2+1 relativistic fermions; 
(ii) lead to the finite density of states P8f , thus trans- 
forming the system to Class (1); (iii) produce zeroes of 
fractional dimension, which means that the exponent in 
the density of states N(E) oc E 2 ~ D is non-integral |2g] 
and thus corresponds to a fractional D of the manifold 
of zeroes; and (iv) lead to localization |3(J. An open 
question is: Can the quantum fluctuations do the same, 
in particular, can they change the effective dimension of 
the zeroes? 



2. Scaling law near zeroes. 

As in the other two classes of fermionic systems with 
the gapless quasiparticles, all low-energy (infrared) prop- 
erties of cuprate superconductors are determined by ze- 
roes. In particular, the density of the fermionic states is 
determined by the dimension of the zeroes: 



N(E) = Y,5(E-E(p))~g 



1-D 



(15) 



Many low-temperature properties of these superconduc- 
tors are obtained from a simple scaling arguments. For 
example, an external magnetic field B has dimension of 
E 2 and thus of T 2 . At finite B, the density of states is 
nonzero even at E = 0. Substituting B ~ E 2 to Eq.jlq) 
one obtains N(0, B) ~ _B( 2 ~ D )/ 2 and the following scal- 
ing law for the heat capacity: 



C(T,B) = B^ 2 - D y 2 Tf(^), 



(16) 



where / is some function with the known asymptotes 
(see fl24|| ). An experimental indication of such scaling 
with D — 1 was reported for YBa2Cu307 in Ref. |2q|. 



III. PROPERTIES OF SYSTEM WITH FERMI 
POINTS. 

A. Relativistic massless chiral fermions. 

Close to the Fermi point p^ ' in the 4D space one can 
expand the propagator in terms of the deviations from 
this Fermi point, p^ — p\ L . If the Fermi point is not 
degenerate, the general form of the propagator is 



^(Pm 



pf) 



(17) 



Here we returned back from the imaginary frequency axis 
to the real energy, so that z = E = —po instead of z = 
ipo; and T a — (1,t). The quasiparticle spectrum E(p is 
given by the poles of the propagator: 



<T(p„ 



pf)(p» 



pl o) ) = 0, g 



flV 



ab n v 
V e a e b 



(18) 



Thus in the vicinity of the Fermi point the massless quasi- 
particles are described by the Lorentzian metric g^ v . It 
is most important that this is the general form of the 
energy spectrum in the vicinity of any Fermi point, even 



if the underlying Fermi system is not Lorentz invariant; 
superfluid 3 He-A is an example. The fermionic spectrum 
necessarily becomes Lorentz invariant near the Fermi 
point, i.e.. If one applies this reasoning to our quan- 
tum vacuum, one may conclude that possibly the ob- 
served Lorentz invariance of the physical laws is not a 
fundamental but a low-energy property of the vacuum 
which results from the existence of the topologically sta- 
ble Fermi points. 



B. Collective modes — electromagnetic and 
gravitational fields. 

Let us consider the collective modes in such a system. 
The effective fields acting on a given particle due to in- 
teractions with other moving particles cannot destroy the 
Fermi point. That is why, under the inhomogeneous per- 
turbation of the fermionic vacuum the general form of 
Eqs.( [17|jl8| ) is preserved. However the perturbations lead 
to a local shift in the position of the Fermi point p^ in 
momentum space and to a local change of the vierbein e^ 
(which in particular includes slopes of the energy spec- 
trum (see Fig. ||). This means that the low- frequency col- 
lective modes in such Fermi liquids are the propagating 
collective oscillations of the positions of the Fermi point 
and of the slopes at the Fermi point. The former is felt 
by the right- or the left-handed quasiparticles as the dy- 
namical gauge (electromagnetic) field, because the main 
effect of the electromagnetic field A^ = (Ao,A) is just 
the dynamical change in the position of zero in the energy 
spectrum: in the simplest case (E—eAo) 2 = c 2 (p — eA) 2 . 

The collective modes related to a local change of the 
vierbein e^ correspond to the dynamical gravitational 
field. The quasiparticles feel the inverse tensor g^ v as the 
metric of the effective space in which they move along the 
geodesic curves 



-d^a. 



These attributes 



ds 2 = g^dx^dx" 



(19) 



Therefore, the collective modes related to the slopes play 
the part of the gravity field (see Fig. g). 

Thus near the Fermi point the quasiparticle is the chi- 
ral massless fermion moving in the effective dynamical 
electromagnetic and gravitational fields. 



" electromagnetic" field A^ — > A u 

of the electromagnetic (A^) and gravitational (g Miy ) fields 

also arise spontaneously as the low-energy phenomena. 

Now let us discuss the dynamics of collective bosonic 
modes, A^ and g^ v '. Since these are the effective fields 
their motion equations do not necessarily obey gauge 
inariance and general covariance. However, in some spe- 
cial cases such symmetries can arise in the low energy 
corner. What are the conditions for that? 

The effective Lagrangian for the collective modes is 
obtained by integrating over the vacuum fluctuations of 
the fermionic field. This principle was used by Sakharov 
and Zeldovich to obtain an effective gravity |3lJ] and ef- 
fective electrodynamics ]32J] - both arising from fluctua- 
tions of the fermionic vacuum. If the main contribution 
to the effective action comes from the vacuum fermions 
whose momenta p are concentrated near the Fermi point, 
i.e. where the fermionic spectrum is linear and thus 
obeys the "Lorentz invariance" , the result of the integra- 
tion is necessarily invariant under gauge transformation, 
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and has a covariant form. The obtained 



effective Lagrangian then gives the Maxwell equations for 
Aft and the Einstein equations for g^, so that the prop- 
agating bosonic collective modes do represent the gauge 
bosons and gravitons. 

Thus two requirements must be fulfilled - (i) the 
fermionic system has a Fermi point and (ii) the main 
physics is concentrated near this Fermi point. In this 
case the system acquires at low energy all the proper- 
ties of the modern quantum field theory: chiral fermions, 
quantum gauge fields, and gravity. All these ingredients 
are actually low-energy (infra-red) phenomena. 

There is another important symmetry obeyed by mass- 
less relativistic Weyl fermions, the conformal invariance 
- the invariance under transformation g^ v — > a(r, t^g^y. 
In the extreme limit when the vacuum fermions are dom- 
inatingly relativistic, the effective action for gravity must 
be conformly invariant. Such gravity, the so-called Weyl 
gravity, is a viable rival to Einstein gravity in modern 
cosmology |33|,[54]]: The Weyl gravity (i) can explain the 
galactic rotation curves without dark matter; (ii) it re- 
produces the Schwarzschild solution at small distances; 
(iii) it can solve the cosmological constant problem, since 
the cosmological constant is forbidden if the conformal 
invariance is strongly obeyed; etc. (see [J35[). 



C. Gauge invariance, general covariance, conformal 
invariance. 



D. He-A: Gauge invariance but no general 
covariance. 



Another important property which results from the 
above equation is that the fermionic propagator in 
Eq.(p7|) is gauge invariant and even obeys the general 
covariance near the Fermi point. For example, the local 
phase transformation of the wave function of the fermion, 
"J — > vfe 16 "^'*) can be compensated by the shift of the 



Let us consider what happens in a practical realization 
of systems with Fermi points in condensed matter - in 
3 He-A. Close to the gap nodes, i.e. at energies £«A, 
where A is the maximal value of the gap in 3 He-A which 
plays the part of the Planck energy, the quasiparticles do 
obey the relativistic equation 



g" v { Pfi -eA fl )(p l/ -eA v ) = 



(20) 



Here e = ± is the "electric charge" and simultaneously 
the chirality of the quasiparticles. Let us consider the 
simplest situation, when the 3 He-A is in its vacuum mani- 
fold, which is characterized by two unit mutually orthog- 
onal vectors e^ and e^ 2 \ and neglect the spin degen- 
eracy of the Fermi point. When we consider the low- 
energy collective modes, these vectors are slowly chang- 
ing in space-time. In this situation the effective metric 
and effective electromagnetic field are given by: 



v%(8 



A = p F \ , A = p F v s ■ 1 , 


(21) 
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g 00 = -l,g 0i = vi, 


(22) 



where v s is the superfluid velocity given by v s i = 
(S/2 m )e( 1 'V«e< 2 ). 

From above equations it follows that the fields, which 
act on the " relativistic" quasiparticles as electromagnetic 
and gravitational fields, have very strange behavior. For 
example, the same texture of the 1- vector is felt by quasi- 
particles as the effective magnetic field B = pf V x 1 
according to Eq.(|2l|) and simultaneously it enters the 
metric according to Eq. (p2|) . Such field certainly can- 
not be described by the Maxwell and Einstein equations 
together. Actually the gravitational and electromagnetic 
variables coincide in 3 He-A only when we consider the 
vacuum manifold: Outside of this manifold they split. 
3 He-A, as any other fermionic system with Fermi point, 
has enough number of collective modes to provide the 
analogs for the independent gravitational and electro- 
magnetic fields. But some of these modes are massive in 
3 He-A. For example the gravitational waves correspond 
to the modes, which are different from the oscillations of 
the 1-vector. As distinct from the photons (orbital waves 
- propagating oscillations of the 1-vector) the gravitons 
are massive p0[ |. 

All these troubles occur because in 3 He-A the main 
contribution to the effective action for bosonic fields come 
from the vacuum fermions at the "Planck" energy scale, 
E ~ A. These fermions are far from the Fermi points 
and their spectrum is nonlinear. That is why in general 
the effective action is not symmetric. 



E. Zero charge effect and Maxwell equations. 

There are, however, exclusions, for example, the action 
for the 1-field contains the term with the logarithmically 
divergent factor ln(A/w) [p0|] . It comes from the zero 
charge effect, experienced by the vacuum of the mass- 
less fermions for whom the 1-field acts as electromagnetic 
field. Due to its logarithmic divergence this term is dom- 
inanting at low frequency to: the lower the frequency the 
larger is the contribution of the vacuum fermions from 



the vicinity of the Fermi point and thus the more sym- 
metric is the Lagrangian for the 1-field. As a result, in the 
very low-energy limit, when the non-logarithmic contri- 
butions can be neglected, the effective Lagrangian for the 
A^ becomes gauge invariant and even obeys the general 
covariance: 



L = 



24tt 2 



In 



A 2 






(23) 



where F^ v is the strength of the effective electromagnetic 
field A^ from Eq.(£l|) and g^ v is the effective gravita- 
tional field from Eq. (p2|) . In this regime the A^ field 
does obey the Maxwell equations in a curved space. 



F. Why 3 He-A is not perfect. 

On the other hand the "Einstein" action for y^ LV is 
highly contaminated by many noncovariant terms, which 
come from the integration over the " nonrelativistic" high 
energy degrees of freedom at " Planck" scale. In this sense 
the 3 He-A, with its given physical parameters, is not a 
perfect model for quantum vacuum. To remove the pol- 
luting noncovariant terms, the integration must be spon- 
taneously cut-off at energies much below the "Planck" 
scale, £« A, for example, due to strong quasiparticle 
relaxation. 

The main reason why 3 He-A is not a good substance, 
is that the Fermi points of the left particles, i.e. at 
p = +Pf1j and the Fermi points of the right particles, 
i.e. at p = — pjrl, are far from each orher in equilibrium. 
The "perfect" condensed matter would be such where all 
the Fermi points are at the origin, at p = 0, as it happens 
in the standard model. However, inspite of the absence of 
general covariance even in the low-energy corner, many 
different properties of the physical vacuum with a Fermi 
point, whose direct observation are still far from realiza- 
tion, can be simulated in 3 He-A. One of them is the chiral 
anomaly 



G. Chiral anomaly. 

The chiral anomaly is the phenomenon which allows 
the nucleation of the fermionic charge from the vacuum 
|p6| ) P7[ . Such nucleation results from the spectral flow 
of the fermionic charge through the Fermi point to high 
energy. Since the flux in the momentum space is con- 
served, it can be equally calculated in the infrared or in 
the ultraviolet limits. In 3 Hc-A it is much easier to use 
the infrared regime, where the fermions obey all the "rel- 
ativistic" symmetries. As a result one obtains the same 
anomaly equation, which has been derived by Adler and 
by Bell and Jackiw for the relativistic systems. The rate 
of production of quasiparticle number n — hr + n^ from 
the vacuum in applied electric and magnetic fields is 



1 
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d^ = ^{el-el)F^F; v 



(24) 



Here ur and ul is the density of the right and left quasi- 
particles; en and cl are their charges; and F* is the dual 
field strength. Note that the above equation is fully "rel- 
ativistic". This equation has been verified in 3 He-A ex- 
periments [Bq,B9[, where the "magnetic" B = ppVxl and 
"electric" E = pf^I fields have been simulated by the 
space and time dependent 1-texture. In particle physics 
the only evidence of axial anomaly is related to the de- 
cay of the neutral pion ir° — ► 27, although the anomaly is 
much used in different cosmological scenaria explaining 
an excess of matter over antimatter in the Universe (see 
review jhj). 



H. Degeneracy of Fermi point as the origin of the 
non-Abelian gauge field. 

In 3 He-A the Fermi point (say, at the north pole) is 
doubly degenerate owing to the ordinary spin a of the 
3 He atom. This means that in cquilibtium the two ze- 
roes, each with the topological invariant N3 — —1, are 
at the same point in momentum space. Let us find out 
what can be the consequences of the Fermi point degen- 
eracy. It is clear that the collective motion can split the 
Fermi points: positions of the two points can oscillate 
separately. Moreover, since the propagator is now the 
4x4 matrix there can be the cross terms. If we neglect 
the degrees of freedom related to the vierbein then the 
collective variables of the system with the doubly degen- 
erate Fermi point enter the fermionic propagator as 



r a e^-eA^-ea a W«) 



(25) 



The new effective field W" acts on the chiral quasipar- 
ticles as SU(2) gauge field. Thus in this effective field 
theory the ordinary spin of the 3 He atoms plays the part 
of the weak isospin @|o[. The "weak" field W* is also 
dynamical and in the leading logarithmic order obeys the 
Maxwell (actually Yang-Mills) equations. 

This implies that the higher symmetry groups of our 
vacuum can in principle arise as a consequence of the 
Fermi point degeneracy. For example, the 4-fold degen- 
eracy of the Fermi point (which implies that there are, 
say, 4 left-handed fermionic species, or 2 left-handed + 
2 right-handed) can produce the SU(4) gauge group. In 
particle physics the collective modes related to the shift 
of the 4- momentum are discussed in terms of the "gen- 
eralized covariant derivative" M3,E3[. In this theory the 
gauge fields, the Higgs fields, and Yukawa interactions, 
all are realized as shifts of positions of the degenerate 
Fermi point, with degeneracy corresponding to different 
quarks and leptons. 

In the Eq.(E5h we did not take into account that dy- 
namically the vierbein can also oscillate differently for 



each of the two elementary Fermi points. As a result the 
number of the collective modes could increase even more. 
This is an interesting problem which must be investigated 
in detail. If the degenerate Fermi point mechanism has 
really some connection to the dynamical origin of the 
non-Abelian gauge fields, we must connect the degener- 
acy of the Fermi point (number of the fermionic species) 
with the symmetry group of the gauge fields. Naive ap- 
proach leads to extremely high symmetry group. That 
is why there should be some factors which can restrict 
the number of the gauge and other bosons. For exam- 
ple there can be some special discrete symmetry between 
the fermions of the degenerate point, which restricts the 
number of massless bosonic collective modes. Another 
source of the reduction of the number of the effective 
field has been found by Chadha and Nielsen pj. They 
considered the massless electrodynamics with different 
metric (vierbein) for the left-handed and right-handed 
fermions. In this model the Lorentz invariance is vio- 
lated. They found that the two metrics converge to a 
single one as the energy is lowered. Thus in the low- 
energy corner the Lorentz invariance becomes better and 
better, and at the same time the number of independent 
bosonic modes decreases. 

Since the connections between QFT in the standard 
model and in 3 He-A has been discussed in the earlier pub- 
lications EG,p4],p9[, we concentrate here on some prob- 
lems related to gravitational analogy. 



IV. BLACK HOLE IN d HE-A FILM. 

A. Gravity by motion of superfluids. Sonic black 
hole. 

As we have seen from Eq.(E2) the gravitational field 
can be simulated in 3 He-A by the motion of the liq- 
uid with the superfluid velocity v s and by the 1-texture 
p4|. In this Section we consider the situation when the 
1- vector is uniform and thus does not produce the effec- 
tive gravitational field, so that the gravitational effects 
come only from the motion of the superfluid vacuu. The 
propagation of fermions in the moving liquid obeys the 
same equation as propagation of relativistic particles in 
the gravitational field. The same happens for propagat- 
ing sound waves in normal fluids pq,Eq| and phonons in 
superfluid 4 He. In the simplest case of the radial mo- 
tion of superfluid He, the effective metric is expressed 
in terms of the radial superfluid velocity v s (r) as 

ds 2 = - (c 2 - v 2 (r)) dt 2 + 2v s (r)drdt + dr 2 + r 2 dfl 2 , 

(26) 

where c now is a speed of sound in 4 He (phonon velocity). 
For 3 He-A fermions the spherically symmetric metric oc- 
curs if 1- field is radial. Then it follows from Eq.(22) 



ds 2 



(c 2 - v 2 (r)) dt 2 + 2v s (r)drdt + —dr 2 + r 2 dn 2 



(27) 



Kinetic energy of superfiow plays the part of the grav- 
itational potential: $ = — v 2 (r)/2. If one chooses the 
velocity field corresponding to the potential of the point 
body of mass M 



v 2 s {r) 



-2$ 



2GM _^r h 

r r 



(28) 



one obtains the Panleve- 

Gullstrand form of Schwarzschild geometry (see e.g. ref. 
Pq ) . Here r^ denotes the position of the event horizon, 
where the velocity reaches the " speed of light" c. If the 
fluid moves towards the origin, the low-energy quasipar- 
tilcles are trapped behind the horizon, since their speed 
c with respect to the fluid is less than the velocity of the 
fluid. 

Such sonic black hole was first suggested by Unruh for 
ordinary liquid J45j . However since all the known normal 
liquids are classical, the most interesting quantum effects 
related to the horizon cannot be simulated in such flow. 
Also the geometry is such that it cannot be realized: in 
such radial flow inward the liquid is accumulated at the 
origin, so that this sonic black hole cannot be stationary. 
In the other scenario a horizon appears in moving soli- 
tons, if the velocity of the soliton exceeds the local "speed 
of light" J44J. This scenario has the same drawback: in 
finite system the motion of the soliton cannot be sup- 
ported for a long time. In a draining bathtub geometry 
suggested in Ref. J47J the fluid motion can be made con- 
stant in time. However the friction of the liquid, which 
moves through the drain, is the main source of dissipa- 
tion. The superfluidity of the liquid does not help much 
in this situation. Horizon does not appear since the "su- 
perluminal" motion with respect to the boundaries of the 
drain tube is not allowed: The superfiow becomes unsta- 
ble and superfluidity collapses (see Eq] ) . Let us suggest a 
scenario, in which this collapse is avoided. The superfluid 
motion becomes quasi-stationary and exhibits the event 
horizon; the life time of the flow state is determined by 
Hawking radiation. 



B. Simulation of 2D black hole 

The stationary black hole can be realized in the follow- 
ing geometry, which is the development of the bathtub 
geometry of Ref. J47j] (see Fig. [|(a)). The superfluid 3 He- 
A film is moving towards the center of the disk, where it 
escapes to the third dimension due to the orifice (hole). 
If the thickness of the film is constant, the flow veloc- 
ity increases towards the center as v s (r) = b/r and at 
r = r] x = b/c reaches the speed of light (now r denotes 



2D "sonic" black hole in 3He-A film 



* Panleve-Gullstrand form of 2D black hole: 

ds 2 = - dt 2 (c 2 -v 2 (r)) + 2v(r) dr dt + dr 2 + r 2 d(|> 2 

* If 3jje-A film is moving to the hole 
v(r) = b / r 

* Horizon is at n, = b / c 
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3He-A fill 
superfluid ^He film 
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(a) 
black hole horizon 



lvl< c 



4He film 




3He-A fdm 



white hole horizon 
(b) 

FIG. 3. Similation of 2D black hole in thin 3 He-A film, (a) 
Draining bathtub geometry, (b) 3 He-A film circulating on the 
top of the 4 He film on torus. 



the radial coordinate in the cylindrical system) . Outside 
the orifice the motion of the liquid is two dimensional 
and the effective metric for the low-energy Bogoliubov 
quasiparticles is 



ds 2 = - (c 2 - v 2 s ) dt 2 + 2v s drdt + dr 2 



c 2 

—^dz 

vi 



(29) 



Here we took into account that the I vector in the film is 
fixed along the normal to the film, so that the "speed 
of light" for quasiparticles propagating along the film 
c ~ 3 cm/sec is much smaller than the Fermi velocity 
vf which corresponds to the "speed of light" for quasi- 
particles propagating along the normal to the film. Note 
that c is much smaller than the speed of sound in 3 He- 
A, that is why the motion of fluid has no effect on the 
density of the liquid. 

The important element of the construction in Fig. is 
that the moving superfluid 3 He-A film is placed on the 
top of the superfluid 4 He film. This is made to avoid the 
interaction of the 3 He-A film with the solid substrate. 
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The superfluid 4 He film effectively screens the interac- 
tion and thus prevents the collapse of the superfluid flow 
of 3 He-A with " super luminal" velocity: Since the inter- 
actions with walls is removed the local observer moving 
with the superfluid velocity cannot detect the relative 
motion with respect to the wall. 

The motion of the superfluid 3 He-A with respect to 
superfluid 4 He film is not dangerous: The superfluid 4 He 
is not excited even if 3 He-A moves with its superluminal 
velocity: c for 3 He-A is much smaller than the Landau 
velocity for radiation of quasiparticles in superfluid 4 He, 
which is about 50 m/sec. In this consideration we ne- 
glected the radiation of surface waves, assuming that the 
thickness of 4 He film is small enough. 

Finally one can close the superflow by introducing the 
toroidal geometry in Fig. |^(b), so that the superfluid 
condensate can circulate. In this case in addition to the 
black hole horizon the white hole horizon appears on the 
path where the superfluid 3 He-A flows out from the ori- 
fice. 

Since the extrinsic mechanism of the friction of 3 Hc-A 
film - the scattering of quasiparticles on the roughness of 
substrate - is abandoned, we can consider now intrinsic 
mechanisms of dissipation. The most interesting one is 
the Hawking radiation related to existence of a horizon. 



C. Vacuum in comoving and rest frames. 

Let us consider the simplest case of the 2D motion 
along the film in the bathtub geometry of Fig. 0(a). 
This can be easily generalized to the motion in the torus 
geometry. 

There are two important reference frames: (i) The 
frame of the observer, who is locally comoving with the 
superfluid vacuum. In this frame the local superfluid ve- 
locity is zero, v s = 0, so that the energy spectrum of the 
Bogoliubov fermions in the place of the observer is (here 
we assume a pure 2D motion along the film) 



En. 



±cp 



(30) 



In this geometry, in which the superflow velocity is con- 
fined in the plane of the film, the speed c coincides with 
the Landau critical velocity of the superfluid vacuum, 
fLandau = rnin(|iJ com (p) \/p). The vacuum as determined 
by the comoving observer is shown on Fig. 0(a) : fermions 
occupy the negative energy levels in the Dirac sea (the 
states with the minus sign in Eq. (|3C|) . It is the counter- 
part of the Minkowski vacuum, which is however deter- 
mined only locally: The comoving frame cannot be de- 
termined globally. Moreover for the comowing observer, 
whose velocity changes with time, the whole velocity field 
v s (r, t) of the superflow is time dependent. This does not 
allow to determine the energy correctly. 



(ii) The energy can be well defined in the laboratory 
frame (the rest frame). In this frame the system is sta- 
tionary, though is not static: The effective metric does 
not depend on time, so that the energy is conserved, but 
this metric contains the mixed component goi — v S i . The 
energy in the rest frame is obtained from the local energy 
in the comoving frame by the Doppler shift: 



-Erest = ±Cp + p-V s (r) 



(31) 



In case of the radial superflow v s (r) = rv s (r) one has 



-E rost = ±cp + p r v s (r) 



(32) 



Figs. [|(b-c) show how the "Minkowski" vacuum of the 
comoving frame is seen by the rest observer (note that 
the velocity is negative, v s (r) < 0). In the absense of 
horizon, or outside the horizon the local vacuum does 
not change: the states which are occupied (empty) in the 
Minkowski vacuum remain occupied (empty) in the rest 
frame vacuum (see Fig. |](b)). In the presense of hori- 
zon behind which the velocity of superflow exceeds the 
Landau critical velocity the situation changes: Behind 
the horizon the vacuum in the rest frame differs from 
that in the comoving frame. Let us for simplicity con- 
sider the states with zero transverse momentum p^ = 
on the branch -E res t = (v s (r) + c)p r in the rest frame. If 
the system is in the Minkowski vacuum state (i.e. in the 
ground state as viewed by comoving observer), quasipar- 
ticles on this branch has reversed distribution in the rest 
frame: the negative energy states are empty, while the 
positive energy states are occupied (see Fig. ||(c)). For 
this branch the particle distribution corresponds to the 
negative temperature T — —0 behind horizon. 

Since the energy in the rest frame is a good quantum 
number, the fermions can tunnel across the horizon from 
the occupied levels to the empty ones with the same en- 
ergy. Thus if the system is initially in the Minkowski 
vacuum in the comoving frame, the tunneling disturbs 
this vacuum state: Pairs of excitations are created: the 
quasiparticle, say, is created outside the horizon while 
its partner - the quasihole - is created inside the hori- 
zon. This simulates the Hawking radiation from the black 
hole. 



D. Hawking radiation 

To estimate the tunneling rate in the semiclassical 
approximation, let us consider the classical trajectories 
p r (r) of particles, say, with positive energy, E rcst > 0, for 
the simplest case when the transverse momentum p^ is 
zero Fig. ||. The branch E rest — (v s (r)—c)p r describes the 
incoming particles with p r < which propagate through 
the horizon to the orifice (or to the singularity at r = 0, 
if the orifice is infinitely small) without any singularity at 



If 



Hawking radiation 




Vacuum in the comoving frame 



E r e S t(r>r h ) 



E resl (r<r h ) 




(V+C)p r 

. (v-c)p r 

(c) f v "C)Pr 
The same vacuum viewed in rest frame: 
(b) outside horizon; (c) behind horizon. 
Tunneling occurs from occupied states behind horizon 
to empty states outside horizon 

FIG. 4. (a) Fermionic vacuum in the comoving frame. 
The states with E corn < are occupied (thick lines). The 
same vacuum viewed in the rest frame (b) outside hori- 
zon and (c) inside horizon. Behind the horizon the branch 
-Erest — (v + c)p r (for p± = 0) has inverse population as seen 
in the rest frame: the states with positive enegry -E rcs t > 
are filled, while the states with -E rcs t < are empty. The tun- 
neling across horizon from the occupied states to the empty 
states with the same energy gives rise to the Hawking radia- 
tion from the horizon. 



the horizon. The classical trajectories of these particles 
are 



Pr(r) 



<0 



(33) 



C — v s (r) 

The energies of these particles viewed by the comov- 
ing observer are also positive: E com (r) — —cp r (r) = 
E Icst {l~{v s {r)/c))- 1 >0. 

Another branch -E rcs t = (v a (r) + c)p r in Fig. ra contains 
two disconnected pieces describing the particle propagat- 
ing from the horizon in two opposite directions: 

-Ercst 



r > r h : p r (r) = 
r <r h : p r (r) = 



c + v s (r) 

-^rcst 

c + v s (r) 



E com {r) = c Pr {r) > (34) 
Scom (r) = c Pr (r) < (35) 



* Panleve-Gullstrand form of black hole: 

ds 2 = - dt 2 (c 2 -v 2 (r)) + 2v(r) dr dt + dr 2 + r 2 d0 2 
kinetic energy of flow = gravitational potential: v 2 (r)/2 = GM/r 

* Particle energy in the rest frame E res t = ± clp[ + p r v(r) 
Particle energy in comoving frame E comov j ng = ± dpi 

* Trajectories of particle with positive energy in the rest frame 
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r=0 



p r = E/(l+v(r)) 
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p r =-E/(l-v(r)) 

Erest = Clprl + PrV(r) > 
Ecomoving = Clprl > U 



particle is in the vacuum state in comoving frame 
but has positive energy E res t in the rest frame 

* Hawking radiation - tunneling of particle 
from the vacuum state within horizon to outgoing state 

W = w e -2S 



S = Im J dr p r = 71 E / lv'(r h )l = E / 2T Hawking 
FIG. 5. Tunneling from Minkowski vacuum within the 
horizon to the outgoing mode. 



The Eq.(p3) describes the outgoing particles - the par- 
ticles propagating from the horizon to the exterior. The 
energy of these particles is positive in both frames, co- 
moving and rest. The Eq.(p5|) describes the propagation 
of particles from the horizon to the orifice (or to the sin- 
gularity). Though for the rest frame observer the energy 
of these particles is positive, these particles, which live 
within the horizon, belong to the Minkowski vacuum in 
the comoving frame. 

The classical trajectory in Eqs.( 34|35| ) is thus disrupted 
at the horizon. There is however a quantum mechani- 
cal transition between the two pieces of the branch: the 
quantum tunneling. The tunneling amplitude can be 
found in semiclassical approximation by shifting the con- 
tour of integration to the complex plane: 



S = Im / dr p r (r) 



exp(-2S) , 

T-Erest 

H{r)\r=r h ' 



(36) 

(37) 



This means that the wave function of any particle in the 
Minkowskii vacuum inside the horizon contains an ex- 
ponentially small part describing the propagation from 
the horizon to infinity. This corresponds to the radiation 
from the Minkowski vacuum in the presence of the event 
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horizon. The exponential dependence of the probability 
on the quasiparticle energy E rest suggests that this radi- 
ation looks as thermal. The corresponding temperature, 
the Hawking temperature, is 



r, 



Hawking 



h\v' s (r)\ r=rh 
2tt 



(38) 



The radiation leads to the quantum friction: the linear 
momentum of the flow decreases with time. This occurs 
continuously until the superfluid Minkowski vacuum be- 
tween the horizons is completely exhausted and the su- 
perfluid state is violated. This leads to the phase slip 
event, after which the number 7V 3 of circulation quanta 
of superfluid velocity trapped by the torus is reduced. 
This process will repeatedly continue until the two hori- 
zons merge. 



E. Discussion. 

The above construction in Fig. ||(b) allows us (at least 
in principle) to obtain the event horizon in the quasi- 
stationary regime, when the main source of nonstationar- 
ity is the dissipation coming from the Hawking radiation. 
As for the practical realization, there are, of course, many 
technical problems to be solved. On the other hand, if the 
black hole analog can in principle exist in condensed mat- 
ter as the quasi-stationary object, its prototype - the real 
black hole - can also exist, at least in principle (though 
it is not so easy to find the scenario of how this object 
can be obtained from the gravitational collapse of matter 

It appears that in some range of the energies of 
Hawking-radiated particles, close to the "Planck" scale, 
where the fermionic spectrum becomes "nonrelativistic", 
the Hawking radiation does not exist any more. The par- 
ticles which are radiated in the fully relativistic case, will 
be now Andreev scattered back to the black hole. Thus 
both partners (particle and hole) of the Hawking radia- 
tion will remain within the horizon, so that the particle 
creation in high gravity field will disturb the quantum 
vacuum inside the horizon without any radiation outside. 
In principle the pair creation inside the horizon can be 
more important for the dissipation in the superfluid 3He 
film under discussion than the Hawking radiation. 

V. ROTATING VACUUM. 
A. Unruh effect. 



so that the matter of the body gets heated to Tjj (see ref- 
erences in [pTJ ) . Linear motion at constant proper accel- 
eration (hyperbolic motion) leads to velocity arbitrarily 
close to the speed of light. On the other hand uniform 
circular motion features constant centripetal acceleration 
while being free of the above mentioned pathology (see 
the latest references in p2-B5|). The latter motion is 
stationary in the rotating frame, which is thus a conve- 
nient frame for study of the radiation and thermalization 
effects for uniformly rotating body. 



B. Zel'dovich-Starobinsky effect. 

Zel'dovich [Q was the first who predicted that the 
rotating body (say, dielectric cylinder) amplifies those 
electromagnetic modes which satisfy the condition 



u> ■ 



LQ < 



(39) 



Here co is the frequency of the mode, L is its azimuthal 
quantum number, and O is the angular velocity of the 
rotating cylinder. This amplification of the incoming ra- 
diation is referred to as superradiance [|57| . The other as- 
pect of this phenomenon is that due to quantum effects, 
the cylinder rotating in quantum vacuum spontaneously 
emits the electromagnetic modes satisfying Eq.(|9|) [|5o[ . 
The same occurs for any rotating body, including the ro- 
tating black hole Kq], if the above condition is satisfied. 

Distinct from the linearly accelerated body, the radia- 
tion by a rotating body does not look thermal. Also, the 
rotating observer does not see the Minkowski vacuum as 
a thermal bath. This means that the matter of the body, 
though excited by interaction with the quantum fluctu- 
ations of the Minkowski vacuum, does not necessarily 
acquire an intrinsic temperature depending only on the 
angular velocity of rotation. Moreover the vacuum of the 
rotating frame is not well defined because of the ergore- 
gion, which exists at the distance r e = c/Q from the axis 
of rotation. 

The problems related to the response of the quantum 
system in its ground state to rotation |5^], such as ra- 
diation by the object rotating in vacuum |5q , |59| , p8| , p7| 
and the vacuum instability caused by the existense of er- 
goregion p0( , etc., can be simulated in superfluids, where 
the superfluid ground state plays the part of the quan- 
tum vacuum. The quantum friction due to spontaneous 
emission of phonons in superfluid 4 He and Bogoliubov 
fermions in superfluid 3 He-B has been discussed in [ pl[ . 
Here we extend an analysis to the radiation of quasipar- 
ticles in 3 He-A. 



The body moving in the vacuum with linear accelera- 
tion a is believed to radiate the thermal spectrum with 
the Unruh temperature Tjj = Ka/2irc |p0|| . The comoving 
observer sees the vacuum as a thermal bath with T = Tjj, 



C. Cylindrical geometry. 

Let us consider a cylinder of radius R rotating with an- 
gular velocity tt in the (infinite) superfluid liquid. In this 
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situation there are again two important reference frames. 
One of them is the laboratory frame. The energy of quasi- 
particles is not well determined in the laboratory frame, 
because, if the rotating body is not the perfect axysym- 
metric cylinder, it produces the time dependent pertur- 
bations of the liquid. However at distances far enough 
from the surface of the cylinder the influence of the ro- 
tating cylinder on the liquid can be neglected. In this ap- 
proximation the superfluid vacuum is in rest with respect 
to the laboratory frame. This means that the superfluid 
velocity v s = in the laboratory frame. The quasiparti- 
cle energy in this frame is that as in the comoving frame 
of previous Section. It is -Ecom = cp if we consider super- 
fluid 4 He, or the Eq. (|30|) , if we consider a pure 2D motion 
of 3 He-A with 1 = z. Such energy spectrum corresponds 
to the effective Minkowski metric of flat space 



ds 2 



-c 2 dt 2 + r 2 d(j) 2 + dr 2 + a 2 dz 2 



(40) 



Here a = 1 for isotropic 4 He and a — c/vf for anisotropic 
3 He-A. 



D. Conical texture with negative angle deficit. 

In the 3D case of 3 He-A there can be another effec- 
tive metric far from the body. It can be caused by the 
normal boundary condition on the 1-vector, which prefers 
the radial orientation of 1-vector. In the case of the radial 
distribution, 1 = f , the effective metric for the quasipar- 
ticles moving outside the cylinder follows from Eq.(E3): 



ds 2 



-c 2 dt 2 + dz 2 + r 2 



c 2 

-^-dr~ 

v F 



(41) 



Such an effective space is conical: The space outside the 
cylinder is flat, but the proper length of the circumfer- 
ence of radius r around the cylinder is not equal to 2irr. 
In the relativistic theories such conical metric can arise 
outside the local cosmic strings, where there is an angle 
deficit. In our case the length of the circumference is 
1~ktvf i 'c, which is much larger than 2irr: This effective 
space exhibits a "negative angle deficit" (for details see 
Ref. f§). 



E. Rotating frame. 

Since the system far outside the rotating cylinder is 
not disturbed by the rotation of the body, the system 
will remain to be in or close to the Minkowski vacuum (or 
vacuum in the conical space in the case of radial 1-vector) 
as viewed in the laboratory frame. However this does not 
hold in the region ajacent to the cylinder, where the su- 
perfluid velocity field is disturbed and time-dependent. 
Also the energy in the laboratory frame is not well de- 
termined because of the time-dependent perturbations. 



The energy is well determined in the frame corotating 
with the cylinder. In the corotating frame the cylinder 
is at rest, and thus the perturbations caused by rotation 
are stationary. The metric in the coorotating frame is 
simplest far outside the rotating body, where the super- 
fluid velocity in the corotating frame is v s = — CI x r. 
Substituting this v s into Eq. (E2) one obtains the interval 
ds 2 = g^udx^dx" , which determines the propagation of 
quasiparticles in the corotating frame: 

ds 2 = -(c 2 - fl 2 r 2 )dt 2 - 2Q.r 2 d<t>dt + r 2 d(f) 2 + dr 2 + a 2 dz 2 

(42) 

Here again a = 1 for phonons in isotropic He and 
a = c/vp for fermionic quasiparticles in anisotropic 3 He- 
A with 1 = z. These metrics correspond to the interval 
in the rotating frame discussed in relativistic theories. 

It is convinicnt to write the quasiparticle spectrum in 
this frame in two different approximations. In classical 
description one has ^corotating = ^com(p) + P v a . In 
the other description, an azimuthal motion of quasipar- 
ticles is quantized in terms of the angular momentum 
L = rp$, while the radial is still treated in the quasiclas- 
sical approximation. In this case the energy spectrum of 
phonons in the corotating frame is 



E, 



corotating 



cp + p-v s = c\ — + p 2 z +p 2 -VlL , (43) 



and the energy spectrum of the Bogoliubov fermions 



-Ecorotating(p) = ±\/ ^ L 2 + C 2 p 2 + V 2 p {p z =p Pf) 2 - CIL . 

(44) 



F. Ergoregion in superfluids. 

The radius r e = c/fi, where goo = in the Eq. (ff2|) , 
marks the position of the ergoplane. In the ergoregion, 
i.e. at r > r e = c/fl, the energy of quasiparticle in 
Eq.(ff4|) can be negative for any L > 1. This means 
that in the ergoregion the Minkowski vacuum is not the 
vacuum for the corotating observer. Situation is similar 
to the case discussed in previous Section. However the 
relevant superfluid velocity is azimuthal now instead of 
radial. As a result there is an ergoplane instead of a 
horizon. 

From the point of view of the corotating observer, the 
Minkowski vacuum is unstable towards the filling of the 
negative energy states in the ergoregion, which means the 
radiation of the quasiparticles from the rotating cylinder. 
However we need a real process, which leads to such radi- 
ation. This radiation can be caused only by the interac- 
tion between the superfluid Minkowski vacuum and the 
rotating object. 
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Let us consider the slow rotations flR <C c. In this 
case the linear velocity of the cylinder at the surface of 
the cylinder, fii?, is much smaller than the Landau criti- 
cal velocity for nucleation of quasiparticles, VLandau — c. 
Thus quasiparticles cannot be nucleated near the surface 
of cylinder. The ergoregion, where \v s \ = Qr > c and 
quasiparticles can be nucleated, is far from the cylin- 
der, r e 3> R. The interaction with the cylinder, which 
produces the matrix element for the radiation, is very 
small. In this situation the most effective mechanism of 
the quasiparticle radiation is the tunneling of quasiparti- 
cles from the liquid ajacent to the surface of the rotating 
body, which plays the part of the rotating detector, to 
the ergoregion. 



G. Rotating detector. 

The simplest rotating detector is the fermionic system, 
which is rigidly connected to the rotating body. In su- 
perfluids the rotating body can be effectively substituted 
by the rigidly rotating cluster of quantized vortices (see 
Fig. ||(a)). Such clusters, that experience the solid-body 
rotation, are experimentally investigated in both phases 
of superfluid 3 He (see e.g. |63|]). The vortex cluster rotat- 
ing in the infinite superfluid liquid represents the vacuum 
in the corotating frame. However this vacuum state is not 
extended to the exterior of the cluster, where one has the 
analogue of Minkowski vacuum. That is why such state 
is not equilibrium, but it is quasistationary metastable 
state which can live for a macroscopically long time. At 
T = the dominating mechanism of the relaxation of 
angular velocity il of the cluster is the process of the 
radiation of quasiparticles. 

Fig. |6|(b) shows the distribution of superfluid velocity 
v s in the laboratory frame. Within the vortex cluster, i.e. 
at r < R, the superfluid velocity, being averaged over the 
vortices, follows the velocity of the solid body rotation 
of the cluster: i.e. < v s >= Ct x r in the laboratory 
frame and thus < v s >= in the frame corotating with 
the cluster. Outside the cluster the superfluid velocity 
decays as Nn/2Trr, where N is number of vortices in the 
cluster and k is superfluid circulation around individual 
vortex. 

Within the cluster one has < v s > = in the coro- 
tating frame, so one can expect that the spectrum of 
quasiparticles, which live within the cluster, is the same 
as the energy spectrum in the frame comoving with the 
superfluid velocity, £ C orotating = E com . In other words 
this spectrum has no — J7L shift of the energy levels, as 
distinct from quasiparticle spectrum far from the body, 
Eqs . ( p3| , p4| ) , which is measured in the same corotating 
frame. This is not exactly true, since the equation 
< v s >= does not imply that v s = locally: it es- 
sentially depends on the position in the vortex lattice. 



Tunneling region 

R < r < r. 



Ergosurface 

r = r e » c/^\ 




Rotating cylinder 
(vortex cluster) 

<V S >=V C |= Q. r in Lab frame: 

vacuum state 

in corotating frame 

(rotating detector) 

Ergoregion: 

r>r e 
vacuum state 

in Lab frame 

(Minkowski vacuum) 



ergosurface: 

V in rotating frame 

reaches 

Landau velocity 

|flr.-v s |=c 



(p r ,L>0) 




Energy spectrum in core Energy spectrum in ergoregion 
(C) 
FIG. 6. (a) Rotating body - rigidly rotating cluster of vor- 
tices. Within the cluster, at r < R, the average superfluid 
velocity equals the velocity v c i of the solid body rotation of 
the cluster. Superfluid state within the cluster plays the part 
of rotating detector: this is the vacuum state in the corotat- 
ing frame. Far outside of the cluster, where velocity of the 
superfluid is zero in the Lab frame, the system is in Minkowski 
vacuum, (b) Distribution of the superfluid velocity in the lab- 
oratory frame. At r — r e f» c/fi the superfluid velocity in the 
corotating frame reaches the "speed of light" - the Landau 
critical velocity v^andau — c. At r > r e there is an ergoregion 
in the corotating frame, where the quasiparticle negative en- 
ergy states are empty, (c) Tunneling of quasiparticles from 
the vacuum state of the rotating detector to the "Minkowski" 
vacuum in the ergoregion. This produces radiation from the 
rotating body (from vortex cluster) and excitation of the ro- 
tating detector (the fermion zero modes in the cores of vor- 
tices) 



The main feature is nevertheless preserved, when one 
considers the spectrum of quasiparticles that live in 
the vortex core. Their energy spectrum in the coro- 
tating state does not depend on rotation velocity f2: 
^corotating = —u>o(p z )Q, where Q is the generalized an- 
gular momentum; u>o(Pz) is the so called minigap, which 
depends on the core structure. This branch crosses zero 
energy level as a function of Q, if one considers this 
quantum number as continuous, and thus represents the 
fermion zero mode (see left part of Fig. 0(c)). The energy 
spectrum of such chiral fcrmions, that live in the vor- 
tex core, has been first calculated by Caroli, de Gennes 
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and Matricon for the Abrikosov vortex in s-wave super- 
conductors fl64J ; on the relation between the number of 
fermion zero modes and the winding number of the vor- 
tex see [pa] . For us it is important that for vortices in 3 He 
the quantum number Q is integer and thus one has the 
states in the detector with zero energy in the corotating 
frame. 



H. Radiation to the ergoregion. 

The radiation of Bogoliubov quasiparticles can be con- 
sidered as the process in which the particle from the zero 
energy state in the detector, -Ecorotating = 0, tunnels to 
the scattering state at the ergoplane, where also its en- 
ergy is £ corota ting = (Fig. |(c)). In the quasiclassical 
approximation the tunneling probability is e~ 2S , where 
at p z = ±pf and flR <§C c: 



S = Im / drp r (r;E = 0) = L 



dr 



1 



1 



Lin- 



R 



(45) 



Thus all the particles with L > are radiated, but the 
radiation probability is smaller for higher L: 



w oc e 



-2S 



2L 



QR 
c 



2L 



luR 
~~cL 



2L 



(46) 



Here u> = flL is energy (frequency) of the radiated quasi- 
particles in the laboratory frame. 

If c is substituted by the speed of light, Eq. (pq) is pro- 
portional to the superradiant amplification of the elec- 
tromagnetic waves by rotating dielectric cylinder derived 
by Zel'dovich ||,||. 

Since each radiated fermion carries the angular mo- 
mentum L, the vortex cluster rotating in superfluid vac- 
uum at T = is loosing its angular momentum and thus 
experiences the quantum rotational friction. The rada- 
tion also leads to excitation of the detector matter. In 
principle the radiation can occur without excitation of 
the detector vacuum, via direct interaction of the par- 
ticles in the Minkowski vacuum with the rotating body 
(cluster) , but the contribution of this process the quan- 
tum friction is smaller [pill . 



I. Discussion. 

The rotational friction experienced by the body ro- 
tating in superfluid vacuum at T = 0, is caused by the 
spontaneous quantum emission of the quasiparticles from 
the rotating object to the "Minkowski" vacuum in the er- 
goregion. The emission is not thermal and depends on 



the details of the interaction of the radiation with the 
rotating body. In the quasiclassical approximation it is 
mainly determined by the tunneling exponent, which can 
be approximately characterized by the effective tempera- 
ture T c ff ~ hil(2/ln(c/Q,R)). The vacuum friction of the 
rotating body can be observed only if the effective tem- 
perature exceeds the temperature of the bulk superfluid, 
T cff > T. For the body rotating with tt = 10 3 rad/s, 
T must be below 10 _8 K. However, high rotation veloc- 
ity can be obtained for clusters. The cluster containg 
two vortices rotate around their center of mass with 
O = k/AttR 2 , where R is the radius of the circular or- 
bit. If the radius R is of order of superfluid coherence 
length, the effective temperature can reach 10 _4 K. 

The process discussed in this Section occurs only if 
there is an ergoplane in the rotating frame. If the super- 
fluid is contained in a finite external cylinder of radius 
i?oxt > R, this process occurs only at high enough rota- 
tion velocity, r e (Q) — c/Q < i?ext, when the ergoplane is 
within the superfluid. On the instability of the ergore- 
gion in quantum vacuum towards emission see also in 
Ref. |0|. 

If r e (fl) > .Rext and ergoregion is not present, then the 
interaction between the coaxial cylinders via the vacuum 
fluctuations becomes the main mechanism for dissipa- 
tion. This causes the dynamic Casimir forces between 
the walls moving laterally (see Review p6[). As in |6(| 
the nonideality of the cylinders, i.e. violation of the rota- 
tional symmetry by the body, is the necessary condition 
for quantum friction. 



VI. DISCUSSION 

In the above examples of the nontrivial space, the ef- 
fective gravitational field acts as a fixed external field. 
The dynamics of this field has not been discussed here. In 
most cases the effective gravity field do not obey the Ein- 
stein equations. This is the main drawback of superfluid 
3 He-A: Since the Fermi points in 3 He-A are too far apart 
from each other, the dynamical equations for the gravita- 
tional field are clumsy. Nevertheless the above textures 
allow us to simulate many phenomena related to quan- 
tum vacuum in the presence of the strong gravity field. 
This is because many properties of the quantum vacuum 
in curved space, which are determined by the geometry, 
do not depend on the dynamical origin of the geometry. 
For example, it is well known that the Hawking radiation 
is a purely kinematic effect and occurs in any geometry, 
if it exhbits an event horizon Eq] . That is why the 3 He-A 
quantum vacuum is a right object for simulation of many 
aspects of physics of vacuum in a curved space. For ex- 
ample what is the effect of the breakdown of the Lorentz 
invariance at higher energy on Hawking radiation. The 
entropy of the black hole can be also investigated using 
the above model, since the microstates within the horizon 
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are well determined and (at least in principle) are com- 
pletely known in the whole energy range including the 
" transPlanckian" region. All these are different aspects 
of the problem of stability of the vacuum in strong grav- 
itational and other fields. The superfluids provide many 
examples of the instability of the superfiuid vacuum, and 
thus allow us to investigate different mechanisms of re- 
laxation of the physical vacuum - the ether. 
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